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Abstract 

We describe a class of evolution systems of linear partial differential equations with the 
Caputo-Dzhrbashyan fractional derivative of order a G (0, 1) in the time variable t and 
the first order derivatives in spatial variables x = (x±, . . . ,x n ), which can be considered 
as a fractional analogue of the class of hyperbolic systems. For such systems, we con- 
struct a fundamental solution of the Cauchy problem having exponential decay outside 
the fractional light cone {(t,x) : |t~ a x| < 1}. 
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1 Introduction 



Evolution equations with fractional time derivatives are among central objects of the modern 
theory of partial differential equations - due both to various physical applications (dynamical 
processes in fractal and viscoelastic media [T6| [T8] ) and to the rich mathematical content of 
this subject; see, for example, the monographs [U EE2] arid references therein. The theory has 
reached sufficient maturity to investigate not only specific classes of equations like fractional 
diffusion and diffusion-wave equations, but to study general systems of fractional-differential 
equations. In particular, Heibig JTU] obtained an L 2 -existence theorem for general fractional 
systems while the author [13] investigated a class of fractional-parabolic systems. 

In this paper we turn to a fractional version of the class of first order hyperbolic systems with 
constant coefficients (for the classical material see [TJ El [2T] ) . A prototype is the fractional 
diffusion-wave equation 



D^w) (t, x) = Au(t, x), t e (0, T],x 



'1.D 



where 1 < < 2, is the Caputo-Dzhrbashyan fractional derivative, that is 



t 



o 

below we explain how to reduce the equation (1.1) to a system of equations of order (3/2 e (0, 1) 
in t and the first order in spatial variables. 

This equation was studied by many authors (0 EJ [HI [15], [191 ED] and others). For the 
Cauchy problem corresponding to the initial conditions 

u(0, x) = Uq(x), u t (0,x) = 0, (1.3) 

a solution is obtained as a convolution with the Green kernel, 

u(t,x)= [ G(t,x-0u (0d£, (1.4) 



where, as it was proved in [7] for n = 1, and in [TjJ] for the general case, 

\G(t,x)\ < cr^ ln {\x\r^ 2 )E{\x\r^ 2 ), 

{1, ifn=l; 
|logz|, if n = 2; 
z- n+2 , ifn>3, 

E(z) = exp (-az^ \ , C, a > 

(here and below the letters C, a will denote various positive constants). 

For a comparison, note that for the wave equation corresponding formally to = 2, with 
the initial conditions (1.3), a counterpart of the kernel G from (1.4) is a distribution supported 
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on the light cone {Ix^ 1 < 1}, if n is even, or its boundary, if n is odd. There is a similar 
behavior, for example, for some symmetric hyperbolic systems of the first order. 

Thus in the fractional case, though the fundamental solution of the Cauchy problem (FSCP) 

is not concentrated on the set j < 1 j, it decays exponentially outside it possessing a kind 

of weak hyperbolicity property. 

In this paper we describe a class of general time-fractional systems of partial differential 
equations with constant coefficients, of the form 



k=i \ X J 



,m, 0<t<T,i6l", 



1.5) 



where < a < 1, 



(*) 



r(i -a) 



s/t' 



V(r) dr - r a ip(0) 



' dx n 



Pjk are first order polynomials with complex coefficients from the derivatives 
(i = v— T), for which, as in the above example, there exists a FSCP with an exponential decay 
away from the set {|x|t~ Q < 1}. This property implies the existence of solutions for the initial 
functions with exponential growth. We call such systems "fractional-hyperbolic" . We follow 
the techniques by Gel'fand and Shilov [8J and Friedman [6], with modifications needed to cover 
the fractional situation. 



2 Auxiliary Results 

2.1. Derivatives of the Mittag-Leffter functions. The Mittag-Leflier function 

OO I 

EJz)=y^—r 4 r (2.1) 

v ; ^ Via +1) V ; 

i=o v ' 

is a fractional calculus counterpart of the exponential function. The function u(t) = E(Xt a ) 
is a solution of the Cauchy problem (O^-u)^) = \u(t), u(0) = 1 (0 < a < 1, A G C); this 
representation of the solution remains valid for matrix-valued solutions in the case where A is 
a matrix. 

The generalized Mittag-Leffler functions 

OO ; 

E aa {z) = V — — r 

^ r(a/ + 7) 

and 

E p ( z ) - V {p)l - 
2^Y{al + 1 )lV 

where 

(p) l = p(p + l)---(p + l-l),l> 1; (p) = l, 
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are also very useful. 

In this paper we will need estimates for derivatives of the function (2.1). It is known (see 
the identities (1.8.23) and (1.9.1) in [TJ) that 



d 



j- k E a {z) = k\E^ l +ak (z 



(2.2) 



and 



KM 



t(p) 



Here i^>\ is the Wright function, 

i^i 



(P,l) 
(7, a) 



CM) 

(7, a) 



T{l + p) z l 
r(a/ + 7 ) l\ 



The asymptotic behavior of 1^1, as z — > 00, was studied in [221 [2]- The results for our case 
are as follows. If | argz| < ^ — e, < e < ™, then 



E p a (z) ~ const ■ exp(z 



z — > 00. 



If I arg(-z)| < (1 - § )tt - e, < £ < (1 - f )vr, then 

#£, 7 (*)~Q(-s)> 

where Q(z) is the sum of residues of the function 

s h-> 2T(-s)r(s + p)/T(as + 7) 
at the points s = — p — v, v = 0, 1, 2, . . .. This means that 



(2.3) 
(2.4) 



a, 



K»| < C\z\^, \z\ > 1, I aig(-*)| < (1 - -)n - s. 



For the remaining region, ^ — e < ± arg z < ^ + e, where £ > is small enough, 

< 7 W~Q(^)+exp(^). 

It follows from (2.2), (2.4), and (2.5) that for any e > 0, 

'Ce Re2l/a , if|argz|<^, 
C e (l + \z\)~ k ~\ if |argz| > <f + e. 



d k 
dz k 



EJz) 



< 



By (2.6), 



«7T 



< C, if I ared > — — . 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



2.2. Eskin's extension of the Paley-Wiener-Schwartz theorem. The Paley- Wiener-Schwartz 
theorem gives an interpretation of a subclass of the class of entire functions of exponential type 
as the set of Fourier transforms of distributions with compact supports. Eskin [5J considered 
the case of entire functions of order greater than 1. 
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Proposition 1 (Eskin). Let Q(s) = Q(s\, . . . , s n ), Si,...,s n G C, be an entire function of 
order p > 1 satisfying the inequality 

\Q(cr + tr)\<C(l + \a\) l e a ^ P , <x,rGM n , (2.9) 

where I is a nonnegative integer. Consider Q(a), er G W 1 , as a distribution from S'(M. n ). Then 
its Fourier transform Q has the form 

N / f) \ 

Q(x) = J2Rk(g-)F k (x) (2.10) 

k=i V x / 

(in the distribution sense) where the orders of differential operators R k do not exceed l + n + 1, 
if this number is even, or I + n + 2 otherwise, F k are continuous functions satisfying, for any 
e > 0, the inequality 

\F k (x)\<C £ e- ( - b -^ P ', k = l,...,N, (2.11) 



I + J, = l,6=4f 
p p p v 



1 / 1 \ P 



ap 



In a further generalization we consider the case where the function Q depends on an ad- 
ditional parameter t G [0, T], is continuous in (s,t), and both the functions Q and O^Q 
(0 < a < 1) have the upper bounds (2.9) with all the parameters independent of t. Then in 
(2.1), the operators R k do not depend on t, the functions F k and H)[ Fk are continuous in (x, t) 
and have the upper bounds as in (2.11) with all the parameters independent of t. This follows 
easily from the proof given in [5] and was noticed (for a = 1) by Friedman ([6], Theorem 3'). 

3 Fractional Hyperbolicity 

3.1. Definition. Let us consider a system (1.5). It is convenient to use its vector form 

(t, x) = P (i-^J U{t, x), < t < T, x G R n , (3.1) 

where U is a function with values in C m , P is a m x m matrix whose elements are first 
order differential operators. In the "dual" representation, P(s) is a matrix whose elements are 
polynomials in the variables Si, . . . ,s n of degree 1. 

Denote by Xi(s), . . . , X m (s) the roots of the characteristic equation 

det(P(s) - XI) = 0. 

Let 

A a (s)= max ReA^ /a (s). (3.2) 

\arg\ k (s)\<aiT/2 

If | argAfc(s)| > an/2 for all k = 1, . . . , m, then we set A a (s) = 0. In any case, A a (s) > 0. 
We call the system (3.1) fractional-hyperbolic, if 

A a (s) < C {\r\ l ' a + log(M + 1)) , s G C n , (3.3) 
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where s = a + ir, a, r e R n . 

3.2. Special classes of systems. Let us consider systems satisfying the condition 

ReAfc(s) < a\r\ +b (k = l,...,m) (3.4) 
where a, b > 0. Then (3.3) is satisfied. Indeed, if | argAfc(s)| < an/2, then 

Re\l /a (s) = |A fe (s)| 1/a cos(-argA fc (s)). 

a 

An elementary investigation shows that cos(^) < (cos y?) 1/,Q \ if \(p\ < an/2. Therefore 

Re \l /a (s) < (Re\ k (s)) l/a < (a\r\ + b) l ' a < C {\r\ l ' a + l) , 
which implies (3.3). 

Note that the condition (3.4) is slightly stronger than the Gel'fand-Shilov hyperbolicity 
condition for differential systems |S]. 

An important class of fractional-hyperbolic systems is that of "symmetric systems" 

dU(t,x) 



(b^U^J (t, x) + J2 A » q + BU(t, x) = (3.5) 
where A v , v — 1, . . . , n are Hermitian matrices, B is an arbitrary matrix. In this case 

P(s) = i \J2s u A u + ib) . 



Let Hk(s), k = l,...,n, be the eigenvalues of the matrix L(s) = ^2 s vA v + iB. Then 

u=l 

ReAfc(s) = Im/ife(s) for each k. By Hirsch's theorem (see Theorem 1.3.1 in Chapter 3 of 

[17J), |/Xfe(s)| < nM where M is the maximum of absolute values of elements of the matrix 
1 n 

—(Us) - Lis)*) = J2 T v A v + K B - B *)- This implies the inequality (3.4). 
2% v= \ 

3.3. Reduction of the diffusion-wave equation to a system. Let us consider the diffusion- 
wave equation (1.1) (where 1 < f3 < 2) with the initial conditions 

u(0, x) = Uq(x), u t (0, x) = Ui(x), (3.6) 

where Uq,U\ are continuous functions. In this section we deal with classical solutions, that is 
we assume that u G C 2 in the spatial variables, u e C 1 jointly in all the variables, there exists, 
for each t > 0, the Riemann-Liouville fractional derivative 

t 

(< ft «) (t, x) = f ^y ^ J(t- r)-'+y T , x) dr, 

o 

and the equalities (1.1) and (3.6) are satisfied pointwise. 
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Denote v(t,x) = (D (/3/2) m) (t,x). On a function ip e C 1 

t 

1 



l)(/3/2) 



*0 (t) 



r(i-f 



, , (t-r)-^V(r)rfr, 



(3.7) 



so that 



pW 2 V) < Ct 1 * -+0, ast^O. 
Therefore v(0, x) = 0, and W@/ 2 'v coincides with the Riemann-Liouville fractional derivative 



(O) 



1 g 
r(i-§)& 



(t-r)-^(r,x) dr, 



In order to calculate D^^v, we write, as in (3.7) 



1-4 



u(t,aO = I 0+ * t u t (t,x) 



where 



r(7) 



(t - r) T_ 1 (p(r) dr, 7 > 0, 



is the Riemann-Liouville fractional integral. Using the identity io+-^o+ = -^o+ +72 (P2], (2.1.30)) 
we find that 



D' > /' / -'„, - —f'h 1 '^, *- —I 2 ' P iu- — fiD> (/,_1) u) 



ft 



x 



5_ 

at 



y ; - [ ,,\ 1 -1- ... u o(a ; ) _ t -p+i 



D o+, t u ) (t, a;) 



d_ 

dt 



u {x) 



so that 

«.«) («. 

Since T^o+j = ^o+.t (PL (2-1-34)), we obtain from the definition (1.2) of that 



r(i-/3) 



h09/2) 



v) (t, x) = (iB^u) (t, x) + 



r(2-/3) 



We have found that the Cauchy problem (1.1), (3.6) is reduced to the inhomogeneous 
problem 



»C8/2) 



t «) (t,x) = Au(t,x) +Mi(x)- r _ 



f(0,x) = 0, «(0, rr) = Uq(x). 



If ui(x) = 0, so that we consider the conditions (1.3), we can reduce the problem further, 
to a system of the form (3.5). Namely, we set 

v°(t,x)= (of /2) «) (t,x), 

■ du(t,x) . 

v J (t,x)= — , j = l,...,n. 

Denoting V(t,x) = (v °(t, x), v 1 (t, x), . . . , v n (t, x)) we obtain the system 

(itf/v) (t ,, ) + £4,^ = 

where 

/ ... -1 ... 0\ 
... ... 



-1 0... 0...0 



\ 0...0 ... 0/ 

(— 1 is in the (j + l)-th place). 

4 Fundamental Solution of the Cauchy Problem 

4.1. Resolvent matrix-function. Looking for a FSCP satisfying (3.1) with the initial condition 
£7(0, x) = S(x), we apply formally the Fourier transform in x getting the Cauchy problem 

(B ( t a) u)(t,s) = P(s)U(t,s), U(p,s) = I, 

where / is the unit matrix. Then 

U(t,s)=E a (t a P(s)). (4.1) 

The function (4.1) is called the resolvent matrix-function corresponding to the system (3.1). 

Below we will denote the inverse Fourier transform (in the sense specified later) of the 
function (4.1) by G(t, x) leaving the notation U (t, x) for a vector- valued solution with a general 
initial vector-function. 

Proposition 2. For all t e (0,T) ; x G R n , 

\\E a (t a P(s))\\ <C{l+t a \s\) m - l e tka{s) . (4.2) 

Proof. Let us obtain an estimate for ||£' a (£ Q; 7- > )|| where V is an arbitrary matrix with 
eigenvalues /ii, . . . ,/i m . Suppose first that all the eigenvalues are different. Then E a {t a V) = 
f(V) where / is the Newton interpolation polynomial 

f(fi) = 61 + b 2 {n - iii) + b 3 (n - ^X/i - /i 2 ) H h b m (n - n x ) •••(// - // m _i), 
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and the coefficients are chosen in such a way that /(//?) = E a (t a fij), j = 1, ... ,771; see, for 
example, [11] regarding functions of matrices. 

It follows from general estimates given in [8] (Chapter II, Section 6.1) that 



\bf.\ < max 

fi£B k 



d^ 1 
kl E a (t a n) 



dfx 



k = 1 m, 



where B k is the smallest convex polygon containing the points . . . , Using (2.7) and (2.8) 
we come to the inequality 

\b k \ < Ct aik ~ 1) e tM , M= max Rea) /a . 

|arg W ]<^ 3 

Therefore 

\\EJTP)\\ = \\f{v)\\ < Ce tM [1 + n\\v\\ + M) + t 2a (\\v\\ + M)(||P|| + N) + • • •] 

< Ce tM (1 + 2t a \\V\\ + ■■■ + {2t a ) {m - 1) \\V\\ m - 1 ) , 

since \fij\ < \\V\\ for all j. Here the constant C does not depend on V and its eigenvalues. For 
continuity reasons, the above inequality holds for any matrix V, not necessarily with distinct 
eigenvalues. 

Now we set V = P(s). Since H^H 2 does not exceed the sum of squared absolute values of all 
the elements, and since polynomials appearing as elements of P(s) have degrees < 1, we find 
that ||P(s)|| < C{1 + \s\). Substituting and recalling that in this case M = A Q (s), we obtain 
the required inequality (4.2). ■ 

If our system (3.1) is fractional-hyperbolic, that is it satisfies (3.3), then it follows from (4.2) 
that 

\\E a (t a P(s)) || < C(l + |a|)«e a *l T l 1/a (4.3) 

with some q > 0. We call the minimal possible nonnegative integer q, for which (4.3) holds, the 
exponent of the system (3.1). In particular, if our system satisfies (3.4), then q < m — 1. 
Note that 

DS Q) £ Q (f*P(s)) = P(s)E a (t a P(s)), 



so that 



Bl a) E a {t a P{s)) < C(l + \a\) q+1 e 



since the matrix P [i^) contains only first order differential operators. 

4.2. The Cauchy problem. Let us consider the system (3.1) satisfying (3.3), with the initial 
condition U(0,x) = U (x), x eW 1 . A vector-function U(t, x), < t < T, x e R™, is called a 
classical solution of the Cauchy problem, if it is continuous in (t, x), as well as its first derivatives 
in x, the fractional integral 

t 



{lfc*U) (t, x) = / (t - r)- a U(r, x) dr 
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has the first derivative in t, continuous in (t, x), and the equation (3.1) and the initial condition 
are satisfied pointwise. 

Though our aim is to construct a classical solution of the Cauchy problem, at the first stage it 
will be defined as a generalized solution. Given the estimate (4.3), the further reasoning is very 
similar to that of [B] (in particular, to the case of correctly posed systems with positive genus 

it); only the exponent — — — , where po is the reduced order, is replaced with . Therefore 

p - fi 1-a 

here and in the proof of the theorem below, we do not repeat the lengthy calculations from [6] 
referring instead to specific sections of Friedman's work. 

Denote by W i (c > 0) the Frechet space of such functions tp G C°°(R) that for any 

d < c, 

\^\z)\ < Cue-V-W^ , j = 0,l,2,..., 
for all z£l. Our basic space $ of test functions is the direct product of n copies of W i . 

1— a 5 

We will need also a larger space $°, the direct product of n copies of W i , < e < c. 

The Fourier transform J 7 maps W i onto the space W^'c of such entire functions ip(s), 
s = a + it, a, t G R, that for any d > -, s G C, 



\s k tp(s)\ < C M exp (a\dT\*) , k = 0,1,2, 



Denote * = J*$, *° = These spaces have an obvious direct product structure. 

It follows from the estimate (4.3) and the description of multipliers in W«>- ([8], Section 
1.2.4) that multiplication by E a (t a P(s)) (0 < t < T) sends * into \1/ , if T is small enough. 
Therefore the above multiplication defines an action between the conjugate spaces (^/°)' — > ty'. 
Thus, if U G ($ )', then % = FU G (* )' C so that £7(i, s) = E a (t a P(s))Uo(s) is a 
generalized solution over Therefore 

L7(t,x) = (G , (t,-)*L7 (-))(x) (4.4) 

is a generalized solution in the sense of the distribution space defined for i < T, if T is small 
enough. 

Theorem. There exists such a positive number 7 that for any initial function Uq possessing 
continuous derivatives D u Uq, \v\ < q + n + 3 (q is the exponent of the system), such that 

\D u U (x)\ < Ce^^, x G R n , (4.5) 
the generalized solution (4-4) ^ s a classical solution of the Cauchy problem, with the estimate 

\U(t,x)\ < Ce 7 ' |x|T ^, 0<t<T,xeM n , (4.6) 
where C, 7' > do not depend on t,x. The fundamental solution G(t, x) has the form 

M / fj \ 

,(t,x) (4.7) 



jxi / a \ 

G(t,x) = J2 R *[dZ)fk 

k=i ^ J 
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(the differential operators Rk of orders < q + n + 3 are understood in the sense of S'(JBL n ) ), ff. 
are continuous functions satisfying the estimates 

|/ fe (t,z)|<CV-^- a *l^\ 71 >0, (4.8) 
with the constants independent oft,x. 

The properties (4. 7)- (4. 8) express the "fractional- hyperbolic" behavior of the class of sys- 
tems considered in this paper. 

Scheme of proof. The representation (4.7) with the estimates (4.8) follow from (4.3) and 
Proposition 1 (including its generalization described in Section 2.2). Note that the inequality 
(2.11) implies the estimate 

|/*(M)l<ae- ( ^~^- £)l * |T ^, 72 >0, 

Taking e < -jT^ 1 ^" we come to (4.8) with 71 = 72/2. 

It is shown in [6], pages 355-356, that if the convolution U{t,x) defined as a distribution 
from $' is in fact a continuous function in (t, x) and has, with all the derivatives appearing 
in (3.1), the upper bounds like that in (4.6) with coefficients 7' < (1 — c^c 1 ^ where c is the 
constant involved in the definition of the space then U is the classical solution. 

Next ([Sj, page 356), in the situation considered here, the convolutions in the distribution 
sense coincide with the classical ones, so that the estimates of convolutions for functions of 
exponential growth and exponential decay (j6], page 362) are applicable, and differential opera- 
tors in expressions like (2.10) can be switched upon smooth convolutors ([6], page 357). Taking 
into account the assumption (4.5) and using the above convolution estimates we obtain (4.6) 
for any given T, if 7 is small enough. ■ 
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